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We have performed temperature dependent tunneling spectroscopy on an impurity-free surface 
area of a LiFeAs single crystal. Our data reveal a highly unusual temperature evolution of su¬ 
perconductivity: at T* = 18 K a partial superconducting gap opens, as is evidenced by subtle, yet 
clear features in the tunneling spectra, i.e. particle-hole symmetric coherence peaks, and a dip-hump 
structure which signals strong-coupling superconductivity. At Tc = 16 K, these features substantiate 
dramatically and become characteristic of full superconductivity. Remarkably, this is accompanied 
by an almost jump-like increase of the gap energy at Tc to about 87% of its low-temperature gap 
value. The energy of the bosonic mode as measured by the distance between the coherence peak 
and the higher-energy dip remains practically constant in the whole temperature regime T < T*. 

The comparison of these findings with established experimental and theoretical results lead us to 
suggest that the bosonic mode is not directly related to incommensurate spin fluctuations that have 
previously been observed in inelastic neutron scattering. 

PACS numbers: 74.55.+V, 74.62.Yb, 74.70.Xa 


The physics of the material LiFeAs seems to differ in 
many aspects from that of canonical iron-based super¬ 
conductors, and accordingly attracts considerable atten¬ 
tion. Unlike the latter, where superconductivity emerges 
from a Fermi surface-nested antiferromagnetic spin den¬ 
sity wave (SDW) state upon doping [1-4], LiFeAs super¬ 
conducts without any doping [5]. The electronic band 
structure of LiFeAs has been much under debate: after 
early de Haas-van Alphen (dHvA) and scanning tunnel¬ 
ing spectroscopy (STS) experiments from which a nested 
Fermi surface had been deduced [6, 7], there exists now a 
comprehensive data set from transport [8], dHvA [9], STS 
[10, 11], and inelastic neutron scattering (INS) [12, 13] ex¬ 
periments which fully support [11, 14] reports of a Fermi 
surface without any nesting from angular resolved pho¬ 
toemission spectroscopy (ARPES) [15-18]. Strong con¬ 
troversy still exists about the nature of superconductiv¬ 
ity: various order parameter symmetries, including sin¬ 
glet 5^ -wave and 5++-wave as well as triplet p-wave 
states have been proposed in both theoretical and exper¬ 
imental work [10, 18-25]. 

In this letter, we report a careful exploration of the 
temperature dependence of the superconducting state of 
LiFeAs as seen in scanning tunneling spectroscopy (STS) 
of a well-defined clean area of the LiFeAs surface. The 
most striking observation is an apparent phase transition 
within the superconducting state with dramatic impact 
on the spectral signatures of superconductivity in the dif¬ 
ferential conductance dl/dV. More specifically, we ob¬ 
serve the onset of faint superconductivity at T* = 18 K 
as evidenced by a particle-hole symmetric depletion and 
coherence peaks seen in the dljdV at zero bias voltage 
and at about ±3-4 mV, respectively. In addition, char¬ 
acteristic strong-coupling dip-hump modifications in the 
dl/dV appear on both polarities at a distance of about 


5-6 meV from the coherence peaks. These features affect 
only a tiny but well resolvable portion of the total dl/dV. 
Upon lowering the temperature (T), however, these fea¬ 
tures become dramatically enhanced at Tc = 16 K and 
rapidly develop into those of a fully gapped supercon¬ 
ducting state upon further cooling. At the transition at 
16 K, the coherence peaks shift to about ±6 mV which is 
already about 87% of the low-T gap value. Remarkably, 
the energy of the coupled bosonic mode stays largely un¬ 
affected from this transition. We discuss this unusual 
non-BCS-like T-evolution of superconductivity in con¬ 
nection with recent ARPES, INS, and theory results. 

In previous experimental works the reported critical 
temperature Tc for stoichiometric LiEeAs scatters be¬ 
tween about 15 K and 18 K [8, 26-34]. Interestingly the 
occurrence of multiple critical temperatures has previ¬ 
ously been reported from successive nuclear magnetic res¬ 
onance (NMR) Knight shift and AC-susceptibility mea¬ 
surements of one LiEeAs single crystal [35]. However, in 
such an experiment, it is very difficult to rule out sample 
inhomogeneity which in principle could cause the men¬ 
tioned probes to respond to different parts of the sample, 
each with a potentially different critical temperature. In 
order to exclude such thinkable complications, we took 
great care in the sample synthesis, its characterization 
and the STS measurements. 

Single crystals of stoichiometric LiEeAs have been 
grown using the self-flux method as described in Ref. 36. 
In order to confirm stoichiometry and homogeneity of 
the sample prior to mounting it into our scanning tun¬ 
neling microscope (STM), we determined [37] the ^^As 
nuclear quadrupole resonance (NQR) frequency and line 
width of the sample as 21.54 MHz and 24 kHz, respec¬ 
tively. While the frequency is consistent with previous 
findings, the measured line width is by a factor of 1.5-3 
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times sharper than that of previously investigated sam¬ 
ples [12, 20, 35, 36] highlighting the extraordinary homo¬ 
geneity of our crystal. Since LiFeAs is highly air sensi¬ 
tive, these steps, and the actual mounting of the sample 
into our home-built STM [38] have been performed in Ar 
atmosphere. Once the sample had been mounted into 
the STM, the sample space was evacuated, and cooled to 
about 5 K. Subsequently, the sample was cleaved in cryo¬ 
genic vacuum just before approaching an electrochem- 
ically etched W-tip for performing the tunneling mea¬ 
surements. Topography measurements have been exe¬ 
cuted in constant-current mode, while for measuring the 
dl/dV we used a lock-in amplifier with a modulation of 
0.4 mV rms at 1.1111 kHz. In the used sign convention 
of data, negative bias voltages probe the occupied states 
of the sample. In all shown spectroscopic data figures, 
an unavoidable offset of -0.6 mV has been corrected. 

During the d//dV-measurements at different T we en¬ 
sured that the tunneling tip always probed the very same 
atoms within a 2 nm x 2 nm area to rule out any possi¬ 
ble influence of sample inhomogeneity. Within this area, 
800 grid spectra (including forward and backward sweep) 
have been recorded and averaged in the range ±35 mV 
with 0.35 mV energy point resolution. Thereby, the sys¬ 
tem was stabilized at -35 mV before and after every single 
forward and backward sweep. At each T-level, we let the 
system sufficiently thermalize in order to have thermal 
drift less than an atom after every set of grid spectra, i.e. 
after about 2 hours. Furthermore, we confirmed an un¬ 
changed tip state before and after each grid spectroscopy 
from comparing corresponding topography scans [37]. 

Fig. la shows a representative topography scan with 
a field of view of 30 nm x 30 nm, measured at 4.8 K. 
The atomically resolved corrugation reveals the position 
of about 3000 surface Li-atoms [10] with a lattice spacing 
of roughly 3.77 A, matching the reported Li-lattice con¬ 
stant of the material [5]. Isolated locations with bright 
contrast represent impurity states of the first layer of the 
material [39]. Beside these, faint structures appear in 
the image, presumably, arising from second layer impu¬ 
rity states as highlighted in Fig. Ic. In order to spec¬ 
troscopically investigate the pristine electronic structure 
of the material we selected a 2 nm x 2 nm area (indi¬ 
cated by a black square area in Fig. la and Fig. lb), 
far away from first and second layer impurity states, and 
performed scanning tunneling spectroscopy as a function 
of T in this area. 

Our results for the averaged differential conductance 
dl/dV within the area of inset of Fig. la as a function 
of applied bias voltage Wias for the different T-levels are 
shown in Fig. 2. Two aspects can be recognized in the 
shown data: i) in the normal state (20 K) the dljdV ex¬ 
hibits a strong asymmetry between occupied and unoc¬ 
cupied states with a hump-like enhancement around zero 
bias, in agreement with previous findings [33, 34]. ii) At 
base temperature of our system (4.8 K) pronounced sig- 



FIG. 1. (Color online) (a) 30 nm x 30 nm area of atomically 
resolved constant current mode topography image of LiFeAs 
{It = 300 pA, Vhias = ±35 mV) measured at T = 4.8 K. 
White arrows indicate the in-plane shortest Fe-Fe directions. 
The atomic corrugation on the surface corresponds to the Li- 
Li (As-As) lattice spacing of 3.77 A. 22 bright impurities from 
the first layer appear within the scan area. Faint signatures 
of impurities presumably of the second layer of the material 
are also visible (green square). Temperature dependent spec¬ 
troscopy has been measured within the black square of 2 nm 
X 2 nm area, (b) Zoom-in into the blue square in Fig. la to 
show atomic contrast in absence of impurities, (c) Zoom-in 
into the green square in Fig. la to show the influence of an 
impurity in the second layer. 


natures of the superconducting state are superimposed 
on the normal state dl/dV background, where the de¬ 
pletion of the dl/dV at zero bias and the appearance 
of coherence peaks at finite bias voltages are the most 
prominent indicators of the superconducting gap. Upon 
increasing T, the data reveal a systematic closing of the 
gap with an apparent Tc of about 16 K (see also Fig. 3a). 
For illustrating the gap closing we plot the dljdV at the 
energy of strongest depletion close to zero bias in the in¬ 
set of Fig. 2. The onset of a strong decrease of dljdV 
below 16 K is clearly visible. However, a close inspec¬ 
tion of the data reveals that the closing of the gap seems 
incomplete even at temperatures T > 16 K. 
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FIG. 2. (Color online) Temperature dependent tunneling 
spectra measured within the black square of Fig. la/b be¬ 
tween 4.8 K and 20 K. The up-arrow indicates the order of 
the curves at Vbias = 0 with increasing temperature. Inset: 
Zero bias differential conductance as a function of tempera¬ 
ture. The horizontal dashed line is a guide to the eye. Vertical 
dashed lines indicate Tc and T*, see text. 


In order to obtain further insight into the data, we 
normalize all tunneling spectra with respect to the nor¬ 
mal state spectrum at 20 K, see Fig. 3b. In this repre¬ 
sentation, all features of the superconducting state be¬ 
come better visible. At small energies within the gap 
(as marked by the position of the coherence peaks) the 
spectra become almost particle-hole symmetric. Remark¬ 
ably, this symmetry does not persist at larger energies. 
In particular, at positive bias a pronounced dip and a 
hump are present in the spectra at all temperatures (see 
Fig. 3b and Fig. 3c). However, at negative bias, the dip 
corresponds to dl/dV values barely below that of the nor¬ 
mal state, whereas the hump structure is even more pro¬ 
nounced. Such dip-hump features are known as the sig¬ 
nature of strong-coupling superconductivity [40], and are 
frequently observed in different Fe-based superconduc¬ 
tors, where these are often interpreted as fingerprints of 
spin fluctuations [33, 41-43]. Surprisingly, the canonical 
spectral signatures of superconductivity, i.e. the coher¬ 
ence peaks, and the dip-hump structure are even clearly 
present in the normalized tunneling data above the afore 
inferred of 16 K. This is unambiguously shown in 
Fig. 3c which focuses on the T-range 16 K to 18 K. Thus, 
in addition to the onset of pronounced superconductivity 
at Tc = 16 K, our LiFeAs sample is in a superconducting 
state that is characterized by faint corresponding spectral 
features already below T* = 18 K. 

After having established this main experimental ob¬ 
servation of this letter, we turn now to thoroughly ana¬ 
lyzing the observed spectral features as a function of T. 
In Fig. 3b we assign the distance between the coherence 


peaks at positive (negative) energy A+ (A“) to the dou¬ 
ble gap value 2 A and the distance between A+ (A“) and 
the dip position at its high energy side to the mode energy 
11+ At base temperature (4.8 K) we find A ^ 6.9 

mV consistent with previous findings [33, 34]. For the 
mode energy we find 11+ ^ |11“| = 5.4T0.1 mV [44]. The 
T-evolution of superconducting coherence peaks, the dip 
positions, and the resulting mode energies are summa¬ 
rized in Fig. 4. The peak positions which provide an es¬ 
timate for the T-dependence of the superconducting gap 
A(T), interestingly, remain almost constant up to almost 
Tc. To be specific, A(Te = 16 K)/A(4.8 K) 0.87. Upon 
increasing T further, A(T) drops abruptly to about 50% 
of its low-T value at 17 K and becomes barely resolvable 
at T* = 18 K. This unusual behavior is in clear contrast 
to any BCS-like weak coupling scenario [45] as has been 
previously suggested by Chi et al. [33]. A very similar T- 
dependence is found for the positions of the dips. These 
remain almost constant at a value of about ±12 mV up 
to Tc = 16 K, and jump-like decrease to about ±9 mV 
at higher T up to T* = 18 K. Remarkably, the mode 
energy exhibits only a small dip at around Tc but stays 
practically uninfluenced at lower and higher T. 

The observation of two transition temperatures Tc and 
T* at an atomically precisely defined microscopic posi¬ 
tion reveals these as an intrinsic property of the mate¬ 
rial. Thereby it corroborates earlier findings where the 
measured critical temperature on the very same sample 
depended on the probing method and offers a reconcilia¬ 
tion of the spread of reported Tc values [ 8 , 26-35]. 

A multiband electronic structure can in principle be 
pictured as a possible source for a complicated super¬ 
conducting state with multiple possible order parameters 
[46]. In LiFeAs, according to the ARPES-derived elec¬ 
tronic band structure [41], the Fermi surface consists of 
quasi two-dimensional hole-like (labelled 7 ) and electron¬ 
like pockets (labelled p) centered around the F- and M- 
points, respectively. Two further hole-like Fermi surface 
pockets (labelled a) are centered around the Z-point. 
The latter are tiny, yet have been reported to possess 
the largest superconducting gap of about 6 meV as com¬ 
pared to about 3.5 to 4 meV at the 7 - and /3-pockets 
[17, 18]. Recent theoretical work which focuses on an¬ 
alyzing possible pairing scenarios on basis of this band 
structure suggests a complicated multigap scenario which 
allows for peculiar T-dependencies [24]. In particular, it 
has been proposed that the a-Fermi surface pockets at 
kz ^ 7T may cause Cooper pairing prior to that of Fermi 
surface pockets at kz ~ 0 , where the a-bands remain be¬ 
low the Fermi level. As a consequence, at temperatures 
just below the onset of superconductivity, the supercon¬ 
ducting state may be very different from that at lower T, 
where all Fermi surface pockets contribute to the super¬ 
conductivity. Considering such a scenario, the observed 
T* = 18 K might be interpreted as the onset of super¬ 
conductivity at kz ~ TT, whereas the critical temperature 
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FIG. 3. (Color online) (a) Waterfall representation of the differential conductance dl/dV for various temperatures. The 
spectrum at 16 K is highlighted in bold, (b) Differential conductance dljdV at various temperatures normalized to that at 
20 K. Black up-arrows indicate the shift of the position of the positive energy dip at towards lower energy upon raising 

the temperature through Tc — 16 K. The down-arrow indicates the coarse position of the negative energy dip at — A“ — . 

(c) Waterfall representation of normalized spectra in b) at 16 K to 18 K. Superconducting coherence peaks and dip positions 
are indicated by solid and dashed vertical lines, respectively. 
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FIG. 4. (Color online) Temperature evolution of super¬ 
conducting coherence peaks (A^, — A“), the dip positions 
(A+ + — A“ — Q“) and the resulting mode energies (Q^ 

and —Q“). 


Tc = 16 K can be related to the onset of full supercon¬ 
ductivity in the complete Fermi surface at all kz. 

The second unusual observation in our data is the 
practically T-independent value of the mode energy 
which implies that the dominant bosonic coupling mode 
remains the same for all T below T*. In view of the 
afore discussed scenario, this suggests that the observed 


mode primarily couples to the a-bands, at all temper¬ 
atures. Previously, the mode has been suggested [33] 
to be connected with incommensurate spin fluctuations 
that have been observed in inelastic neutron scattering 
[12, 13]. However, the latter have been shown to be un¬ 
related to particle-hole excitation of the a-bands but in¬ 
stead emerge from those between the 7 - and the /5-bands 
[13, 14]. Thus, the observed bosonic mode is likely to 
stem from another, yet to be discovered excitation in the 
system. 

In conclusion, we have presented temperature depen¬ 
dent tunneling spectroscopy on the bare surface of sto¬ 
ichiometric LiFeAs with stable atomic registry. We ob¬ 
serve a highly unusual temperature evolution of the su¬ 
perconducting state with the onset of faint superconduc¬ 
tivity at T* = 18 K, which develops into full supercon¬ 
ductivity only at Tc = 16 K. In the whole superconduct¬ 
ing temperature regime T < T* we observe evidence for 
a strong coupling bosonic mode that remains constant in 
energy except a slight reduction at Tc. The comparison of 
our findings with results from inelastic neutron scatter¬ 
ing, experimental band structure data from ARPES, and 
theory let us to suggest that the bosonic mode might be 
unrelated to observed incommensurate spin fluctuations. 
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